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Abstract 


In this article we will only touch on a few tiny parts of the field of linear Dio- 
phantine equations. Some of the tools introduced, however, will be useful in many 
other parts of the subject. 


1 Whole Numbers 


In number theory, we are usually concerned with the properties of the integers, or 
whole numbers: Z = {..., —3, —2, —1,0,1,2,3,...}. Let us begin with a very simple 
problem that should be familiar to anyone who has studied elementary algebra. 

e Suppose that dolls sell for 7 dollars each, and toy train sets sell for 18 dollars. A store 
sells 25 total dolls and train sets, and the total amount received is 208 dollars. How 
many of each were sold? 

The standard solution is straight-forward: Let x be the number of dolls and y be the 
number of train sets. Then we have two equations and two unknowns: 


arty = 25 
7x+18y = 208 
The equations above can be solved in many ways, but perhaps the easiest is to note that 


the first one can be converted to: x = 25 — y and then that value of x is substituted into 
the other equation and solved: 


7(25—y)+18y = 208 
175—T7y+18y = 208 
—Ty+18y = 208-175 
lly = 33 
y = 3 


Then if we substitute y = 3 into either of the original equations, we obtain x = 22, 
and it is easy to check that those values satisfy the conditions in the original problem. 
Now let’s look at a more interesting problem: 

e Suppose that dolls sell for 7 dollars each, and toy train sets sell for 18 dollars. A store 
sells only dolls and train sets, and the total amount received is 208 dollars. How many 
of each were sold? 


This time there is only one equation: 7x + 18y = 208. We probably learned in algebra 
class that you need as many equations as unknowns to solve problems like this, so at 
first it seems hopeless, but there is one additional key piece of information: the number 
of dolls and the number of train sets must be non-negative whole numbers. With that in 
mind, let’s see what we can do, ignoring for the moment the fact that we already have 
a solution, namely: x = 22 and y = 3. 


Beginning from our original equation, we can do this: 


7x = 208— 18y 
208 — 18y 5 — 4y 
= —— = 29 —-2y4+ . 
- 7 ei Ae 

At first, this seems like we haven’t made any progress, but notice that since y has to be 
an integer, the value of (29 — 2y) is a whole number, so also the fraction (5 — 4y)/7 
also has to be a whole number. Let’s call it a, and then do some arithmetic with the 
resulting equation: 


5 — 4y 
an: 
7a = 5—4y 
4y = 5-—Ta 
is BSG 2 eee 
4 4 


As before, we still don’t have a solution, but things look better in a sense: the denomi- 
nator in the fraction is smaller: it is now 4 instead of 7. As before, since we know that 
a is a whole number, we know that 1 — a is a whole number, and therefore (1 — 3a) /4 
must be a whole number which we will call 6, and we’ll again repeat the same sort of 
algebraic manipulations: 


a 1-—3a 
4 
4b = 1-3a 
3a = 1-46 
1 — 4b ee 
a= —=- —-.. 
3 3 


The same reasoning as before tells us that (1 — b)/3 must be a whole number, so we 
call it c: 


£. 26 
SS 
3c = 1-5 
b = 1-38c. 


Now we’ve actually made some progress. No matter what integer value c takes, b will 
be an integer! Let’s substitute that value of b back into the previous equation: 


(adh, eBay. 3 215. 
3 3 3 ire 


a 


Now substitute this value of a to obtain y: 


6-71 5—7(-1+4c) _ 12—28c _ 


Finally, we can substitute this value of y into the original equation to obtain 2: 


208 — 1 208 — 18(3 — 154+ 12 
= ae 8y _ 08 se Te) _ a Oe 90418 


Our solution looks a little different from what we obtained in the first problem, but here 
it is: 


22 + 18¢ 
= 3-7c 


8 
I 


If c = 0 we obtain the same solution we did previously: x = 22 and y = 3, but note 
that any integer value of c will yield another solution. We can see that positive values 
of c will yield negative values of y, but if c = —1, we obtain another solution: « = 4 
and y = 10. It’s easy to check that both (2, y) pairs are valid solutions to the original 
problem. If cis smaller, —2 or less, then the x values become negative, so there are no 
additional solutions. 

When an equation of this sort is solvable by this method, there is no limit to the number 
of steps that need to be taken to obtain the solution. In the example above, we needed 
to introduce integers a, b and c, but other equations might require more or fewer of 
these intermediate values. 


2 Linear Diophantine Equations 


What we have just solved is known as a Diophantine equation — an equation whose 
roots are required to be integers. Probably the most famous Diophantine equation is 
the one representing Fermat’s last theorem, finally proved hundreds of years after it 
was proposed by Andrew Wiles: 


If n > 2, there are no non-trivial! solutions in integers to the equation: 
gh + y” _ Pan 


There are many, many forms of Diophantine equations, but equations of the sort that 
we just solved are called “linear Diophantine equations”: all the coefficients of the 
variables are integers. 

Let’s look a little more closely at the equation we just solved: 7x + 18y = 208. If the 
only requirement were that the roots be integers (not necessarily non-negative integers), 
then our solution: x = 22 + 18c and y = 3 — 7c represent an infinte set of solutions, 
where every different integer value of c generates another solution. 


'The “trivial” solutions are those where x or y is zero. 


A more geometric view of the problem is this: If we were to graph the equation 7x + 
18y = 208, the solutions are places where the graph passes through points that have 
integer coordinates. In Figure 1 a portion of that line is plotted, and the points where 
the graph has integer coordinates are indicated and labeled. 


(-32,24) 
aes 


(410) 


(22,3) 


_(40,-4) 


Figure 1: Graph of 7x + 18y = 208 


Notice that all the points with integer coordinates are evenly spaced along the line. 
In fact, if we begin at any point and add 18 to the x-coordinate and at the same time 
subtract 7 from the y-coordinate, we arrive at another point on the graph with integer 
coordinates. A quick examination of the original equation should make it obvious why 
this is the case. The equation is: 


7x + 18y = 208. 


If we add 18 to the x value, we increase the left side by 7 - 18. If we subtract 7 from 
the y value, we decrease the left side by the same amount: 18 - 7. The net effect is to 
leave the left side unchanged. 

Notice that this line has a negative slope and happens to cut through the first quadrant 
(quadrant I) and intersect some points with integer coordinates there. This may or may 
not be the case for the graphs of other linear Diophantine equations. Lines with positive 
slopes can have an infinite number of solutions where both are positive, and there are 
equations where there are none. It’s easy to construct such equations with whatever 
characteristics you wish. 


Does every equation of the form: 
ax + by =c, 


where a, b and c are integers have a solution (x, y), where x and y are also integers? 
The answer is no. For example, what if a and b are even and c is odd? The left side 
must be even, and if the right side is odd, there is no possibility of a solution with 
integer values. Similarly, if a@ and b are both multiples of 3 and c is not, the left side 
will be a multiple of 3 and the right side is not, so again, there are no possible integer 
solutions. 

In fact, if the greatest common divisor (GCD) of a and 6b does not divide c, then there 
are no integer solutions. The amazing thing, however, is that if the GCD of a and b also 


divides c, then there are an infinite number of integer solutions, and we will see why 
that is the case later on. 
Note also that another observation we made about our particular problem will also 
apply to a general linear Diophantine equation; namely, that if (a, y) is an integer 
solution to: 

ax +by=c 


then so will be (x + bk, y — ak) where k is any integer. If we substitute x + bk for x 
and y — ak for y, we obtain: 


a(x + bk) +b(y—ak) = c 
ax+abk+by—abk = c 
ax+by = 6, 


so if (x, y) is a solution, then so also is (a + bk, y — ak). 


3 Reducing Fractions 


Now we will jump to a problem that at first glance is totally unrelated to linear Dio- 
phantine equations, but will turn out to be very similar. Let’s begin with a very ugly 
problem: 


Reduce the following fraction to lowest terms: 


179703941 
215237573: 


Recall that to reduce a fraction to lowest terms, you search for numbers that are com- 
mon factors of the numerator and denominator and if such numbers exist, the fraction 
can be simplified by dividing both numerator and denominator by that common num- 
ber. For example, in the fraction 48/66 both the numerator and denominator have a 


factor of 6, so we have: 
48 6-8 8 


66 6-11 11 
It’s easy to check that the 8 and 11 in the fraction 8/11 have no other integer common 
factors, so 8/11 is reduced to its simplest possible form. 
With a computer, or with a heck of a lot of effort by hand or even with a calculator, 
we can search for common factors of the numerator and denominator of our original 
problem, but that could require a great amount of effort. In fact, you will need to test 
hundreds of possible factors before you arrive at the following: 


179703941 — 185071 -971 971 


215237573 = 185071-1163 1163’ 


and even then, a bit more work must be done to assure that 971 and 1163 have no 
common factors. 


What we would like to do is have a method to find the GCD (greatest common divisor) 
of the numerator and denominator, and then eliminate that from the numerator and 
denominator, yielding the fraction reduced to simplest form. 


Let’s try the following approach, which, except for some slightly ugly arithmetic, 
should seem quite familiar: 


215237573 = 179703941 - 1 + 35533632. 


Notice that if some number does divide 215237573 and 179703941 then it must also 
divide 35533632. That means that: 


GCD(215237573, 179703941) = GCD(179703941, 35533632). 


Thus we have reduced the size of the numbers in our problem. 


We can continue in the same way: 


179703941 = 35533632 -5 + 2035781 
35533632 = 2035781 -17 + 925355 
2035781 = 925355-2 + 185071 
925355 = 185071-5+0. 


In every case, any number that divides the two leftmost numbers in the equations above 
must also divide the rightmost, and in the final line, we observe that 185071 must be 
the GCD of those numbers. Once we know that 185071 divides the numerator and 
denominator and in fact is the largest number to do so, we can divide the numerator 
and denominator of our original fraction by that number to obtain the form 971/1163 
as the simplest form. 


This method to obtain the GCD of any two numbers is known as Euclid’s algorithm. 


4 Euclid’s Algorithm and Diophantine Equations 


Now let’s use the Euclidean algorithm on two of the numbers from the original Dio- 
phantine equation we solved in Section 1: 7x + 18y = 208. 


18 = 7-2+4 
7 = 4-143 

= 3*1l+1 
3 = 1-340 


In this example, the final number is 1, so the GCD of 18 and 7 is 1 (in other words, 
18 and 7 are relatively prime), but the interesting thing to note is that the numbers in 
the GCD calculation: 18,7, 4,3, 1 are the same numbers that we got as denominators 
and as the coefficients of the variables in the numerators in the fractions when we were 
solving the Diophantine equation 7x + 18y = 208. The only oddball numbers were 


the constants in the numerators, and that’s not surprising: we never used the number 
208 when we were using the Euclidean algorithm to find the GCD of 18 and 7. If you 
check the arithmetic calculations that are being done in each case, it will be clear why 
the numbers generated in both examples must be the same. 

Suppose that the original Diophantine equation had had a 1 instead of the 208. To make 
sure you understand the technique we used to solve our Diophantine equation it would 
be a good exercise to solve the following equation by yourself before reading on: 


7x + 18y =1 


Here’s the solution (but just the equations: the textual arguments are omitted): 


7a+18y = 1 
1-1 1-4 
—— PN a y 
7 
1—Ay 
a = 
7 
1— 7a jp bese 
= =—-a 
: 4 
—_ 1-— 3a 
4 
1 — 4b one 
a => CO —— 
3 3 
1-6 
c= — 
3 
b = 1-3c 


The nice thing about the 1 in place of the 208 is that it remains constant throughout 
the calculation, whereas the 208 was reduced as various of the denominators divided 
it evenly. In this calculation, all the other coefficients are the same as the numbers 
generated in the straight-forward calculation of the GCD of 7 and 18. 

To complete the solution, we need to back-substitute the b = 1 — 3c and after a few 
steps we obtain: « = —5 + 18c and y = 2 — 7c, where c is an arbitrary integer. 
(Obviously this equation will have no solutions where both x and y are positive.) 
Thus when you do a GCD calculation of a and b, and that GCD turns out to be 1, you’ve 
done a lot of the work toward solving the Diophantine equation 


ax + by = 1. 


So if we can do the Euclidean algorithm, we can find with almost no effort other than 
a little arithmetic the coefficients we need to solve a linear Diophantine equation of the 
form ax + by = 1. Of course we’d like to be able to solve equations where the 1 is 
replaced by an arbitrary number c, but that is actually not too difficult. 


As as example, let’s find solutions for 7z + 18y = 208 assuming that we’ve solved 
7x + 18y = 1. The solutions for the latter equation are x = —5 + 18c and y = 2 — 7c, 


where c is an arbitrary integer. An easy solution is simply to set c = 0 and obtain 
x = —5and y = 2 asa particular solution. But if we multiply xz and y by 208, then the 
left side will be increased by a factor of 208 so if we increase the right side by the same 
factor, we'll have an (x, y) pair that satisfies our original equation 7x + 18y = 208. 
Thus a solution is this: « = —5 - 208 = —1040 and y = 2- 208 = 416. It’s easy to 
plug these numbers in to check that they are valid. 

But we also noticed that adding any multiple of 18 to x while at the same time adding 
that same multiple of —7 to y will yield the other solutions, so the general solution to 
our original problem is: 2 = —1040+ 18k and y = 416 — 7k. If k = 58, for example, 
this yields the solution x = 4 and y = 10. 

We have seen that if we have any solution to one of these linear Diophantine equations, 
we can obtain all the others by adding constant multiples of the opposite coefficients to 
the given solution, all we really need is one solution. 

In the previous examples, once we got to the point where we had b = 1 — 3c, we back- 
substituted and carefully kept track of the coefficient of c in the calculations. But since 
any solution will generate all the others, why not let c = 0? Then we just need to track 
a single number. 


5 Putting It All Together 


Let’s use the techniques above, but in their most simplified form, to solve another 
Diophantine equation. Here’s the problem: 


e In a pet shop, rats cost 5 dollars, guppies cost 3 dollars and crickets cost 10 cents. 
One hundred animals are sold, and the total receipts are 100 dollars. How many rats, 
guppies and crickets were sold? 


If r, g and c represent the number of rats, guppies and crickets, respectively, we’ve got 
two equations (but three unknowns): 


r+g+e = 100 
5r+3g+.le = 100 
To turn the problem into a purely integer problem, multiply the second equation by 10: 
r+gte = 100 
50r+30g+e = 1000 


If we subtract the first equation from the second we obtain the familiar looking linear 
Diophantine equation in two variables: 


49r + 29g = 900. 


Luckily, the GCD for 49 and 29 is 1 which divides 900 so there will be solutions 
(although possibly not solutions where all the values are non-negative. 
(This problem is probably much easier to solve using “guess and check” techniques: 


we know that the number of crickets must be a multiple of 10, so you could just try 0, 
10, 20, ..., 100 of them.) 


Let’s find the GCD of 49 and 29, using the Euclidean algorithm: 


49 = 29-1+20 
29 = 20-1+9 
20 = 9-242 
9 = 2-441 
2 = 1-2+0 


Since we’ve used the variable c for the number of crickets, and g for the number of 
guppies, let’s use variables 7, 7, k, et cetera for the integer values of fractions that we 
get as we step through the solution of the Diophantine equation. We will a particular 
solution to the simpler equation: 


49r + 29g =1 


and then multiply both r and g by 900 for a solution to the original problem. 

Without doing any calculations, but just reading the values obtained from the Euclidean 
algorithm used to calculate the GCD of 49 and 29 and the variables 2, 7, & and / as 
the integer values of fractions, we can just write down the relationships among them 
(expressed both ways). For the purposes of finding a solution to the equation, only the 
expressions on the right are important. Note how the numbers in the fractions (other 
than the 1, of course) are exactly the same as the numbers on the left in the execution 
of the Euclidean algorithm to find the GCD of 49 and 29 above. 


_ 1-299 _ 1-49r 
= gi Hare oe 
_1-20r 1-294 
eae an 
a ee 2. A 0F 
ot 4 o> $< ——|— 
20 9 
1-2) 1—9k 
(eee god! 
9 2 
1-—k 
2 


From the last line, clearly & = 21 for arbitrary integers /, so just set ] = 0 to obtain a 
particular solution. If / = 0, then& = 1. If k =1 then 7 = —4. If 7 = —4, thenz = 9. 
If i = 9 thenr = —13. And if r = —13, then g = 22. 


It’s easy to check our work; namely, that r = —13 and g = 22 is a solution to: 
49r + 29g = 1. 
To obtain a particular solution to the original equation, multiply by 900: r = —11700 


and g = 19800. From previous considerations, we know that the general solution to 


the original equation will be: 


r = —11700+ 29k 
g = 19800 — 49k 


We’re looking for non-negative values of r and g, so divide 29 into 11700 and we find 
that if k = 404 we obtain the values r = 16 and g = 4. If k = 403 or k = 405, either 
r or g is negative. 

Since there are 100 total animals, c = 80, and it’s easy to check that r = 16, g = 4 and 
c = 80 solves the problem. 
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